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Abstract 

We consider a generic mechanism via which thermal relic WIMP dark matter may be 
decoupled from the Standard Model, namely through a combination of WIMP annihilation to 
metastable mediators with subsequent delayed decay to Standard Model states. We illustrate 
this with explicit examples of WIMPs connected to the Standard Model by metastable bosons 
or fermions. In all models, provided the WIMP mass is greater than that of the mediator, it 
can be secluded from the Standard Model with an extremely small elastic scattering cross- 
section on nuclei and rate for direct collider production. In contrast, indirect signatures from 
WIMP annihilation are consistent with a weak scale cross-section and provide potentially 
observable 7-ray signals. We also point out that 7-ray constraints and flavor physics impose 
severe restrictions on MeV-scale variants of secluded models, and identify limited classes 
that pass all the observational constraints. 
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1 Introduction 



The overwhelming astrophysical and cosmological evidence for dark matter has in recent 
years led to a dramatic expansion in experimental programs that aim to detect observa- 
tional signatures of its non-gravitational interactions pp. These probes range from direct 
production at colliders to the recoil of galactic dark matter on nuclei in underground detec- 
tors, and indirect detection of annihilation products, primarily 7-rays. A driving paradigm 
in developing these probes is that of WIMP (weakly interacting massive particle) dark mat- 
ter, which represents a simple and attractive candidate through the fact that a thermal relic 
with weak scale mass and annihilation cross-section into Standard Model (SM) states natu- 
rally provides roughly the correct cosmological abundance [2]. This weak-scale annihilation 
cross-section, when reversed, naturally suggests a weak-scale production cross-section at col- 
liders, and when viewed in the t-channel implies an elastic scattering cross-section on nuclei 
which is within reach of purpose-built underground detectors. In recent years, experiments 
of the latter type have reached an impressive level of sensitivity [21111, and significant future 
progress in this direction is anticipated. 

Along with direct scattering on nuclei, and the neutrino signal from the annihilation of 
WIMPs inside the solar core, which is sensitive to the WIMP trapping rate (and again to the 
elastic scattering cross section), the indirect detection of WIMPs via the products of their 
annihilation in the center of the galaxy is a distinct possibility [5]. This annihilation signal 
can be very model dependent, due to uncertainties in the halo profile, varying branching 
ratios, and the presence or absence of detectable monoenergetic photons. Thus, prior to 
specifying a particular model, it is impossible to say which strategy, direct or indirect detec- 
tion, will provide a more sensitive probe of WIMPs. However, within the prevailing WIMP 
paradigm as exemplified by the lightest superpartner (LSP) in the MSSM for example, the 
stringent constraints on the elastic scattering cross-section often impose significant limits on 
the available indirect signal from annihilation. 

In this paper, we will revisit this aspect of WIMP physics, and question the commonplace 
assumption of a close link between the cross-sections relevant for direct and indirect detec- 
tion. More precisely, we consider what constraints the picobarn annihilation cross-section 
required of a thermal relic WIMP can actually impose on its interactions with normal mat- 
ter, e.g. production at colliders or scattering off nuclei. We observe that in relatively simple 
models, the latter interactions can be highly suppressed and thus in many cases the indirect 
annihilation signature will be the most important probe of the non-gravitational interactions 
of dark matter. 

Generically, any WIMP dark matter model can be conveniently decomposed into three 
sectors, the SM, the WIMP itself, and the fields which mediate the WIMPs interations with 
the SM, 



and in many models the mediator states are in fact part of the SM, e.g. the electroweak gauge 
bosons or the Higgs. In this paper, we point out a simple and generic mechanism that allows 
the WIMP, while still a thermal relic, to be secluded from the SM, dramatically reducing its 
couplings to SM states, and consequently suppressing the collider and direct detection rates 
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Figure 1: The secluded WIMP dark matter scenario. 



by many orders of magnitude. The mechanism relies on a metastable mediator that couples 
the SM to the secluded WIMP sector, with a mass less than that of the WIMP, see Fig. [TJ 
In this kinematic regime, direct annihilation into a pair of mediators is always possible. Due 
to their coupling to the Standard Model, these particles are unstable but the constraint 
on their lifetime is very weak, and in particular a lifetime under one second is sufficient to 
guarantee their decay before the beginning of primordial nucleosynthesis (BBN), rendering 
them completely harmless. Secluded WIMP models may therefore be impossible to detect 
using colliders or direct searches, but the indirect signatures, through e.g. annihilation in 
the Galactic center, can be as pronounced as in any WIMP scenario. 

To illustrate this mechanism, in Section 2 we will construct several models with fermionic, 
scalar and vector particles as mediators. We show that if mwmp < ^mediator, the parameter 
space of such models is highly constrained, as the coupling of the mediator to the SM must 
necessarily be sizable to ensure the required annihilation cross-section. Yet if the reverse is 
true, mwiMP > ^mediator, there are no strict requirements on the size of the mixing except 
for the lifetime of the mediator state, which in some instances can be satisfied for (mixing) 2 
of the mediator with the SM as low as 10~ 23 . 

An interesting limiting regime of the secluded scenario arises when the mediator (and 
the WIMP) are both very light. This ties in with models of MeV-scale dark matter that 
provide a tantalizing yet speculative link [6] to the unexpectedly strong 511 keV line observed 
from the galactic center |7j. We show in Section 3 that some secluded WIMPs may have 
advantages over existing models of this type, and discuss their observational signatures at 
low energies, including missing energy signals in meson decays. 



2 Models of secluded WIMPs 

For some time the discussion of WIMPs has centered around supersymmetric models where 
the LSP is stable, provided that i?-parity is unbroken. Supersymmetry is largely motivated 
by a well-known combination of theoretical arguments unrelated to dark matter, and the 
possible existence of a stable or long-lived WIMP-LSP may provide an interesting bridge 
to cosmology. On the other hand, to date there are no experimental indications of super- 
symmetry, while there is ample evidence for the existence of dark matter. Therefore, an 
alternative approach to the particle physics of dark matter that is certainly logical and jus- 
tifiable, consists of studying generic classes of WIMP models, among which the minimal 
choices are obviously well motivated. Over the years, WIMP models with Higgs and/or 
singlet mediation have been studied extensively [8l [10] . More recently, models with exotic 



electroweak matter were also considered in some detail 12 , as well as models with addi- 
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tional gauge groups [IB] . More generally, going beyond the WIMP framework also allows for 
freedom in tuning the coupling to the SM, as in scenarios with sterile neutrinos [llj or dark 
matter populated by late decays [2] . In this paper we adhere to a rather minimalist WIMP 
framework, which is well suited to demonstrating our main point. 

2.1 U(iy mediator 

We can construct a simple secluded model, starting from the Lagrangian for a Dirac WIMP, 
whose interaction with the SM is mediated by an additional U(l)' gauge group: 

Ik, 

AviMP+mcdiator = ~J V ^u ~ ^V^B^ - \D ^ - U (00*) + 1p(lD ^ - m^jlj) . (2) 

In this Lagrangian, ip and denote respectively the Dirac WIMP and the U(l)' vector 
boson mediator, with field strength V^ u and covariant derivative = <9 M + ie'V^. To avoid 
confusion, we denote the strength of the U(l)' coupling constant as e'. The U(l)' vector 
bosons V M couple to the SM hypercharge gauge bosons via kinetic mixing. The additional 
scalar higgses U(l)' at or near the weak scale. Note that the SM and WIMP sectors are 
coupled only via the mediator at the renormalizable level, the SM is neutral under U(l)', 
and the WIMP sector is a singlet under the SM gauge group. Assuming the scalar potential 
U has a minimum away from zero, the Lagrangian below the U(l)' breaking scale takes the 
form, 

AviMP+mediator = ~\ V ^v + ^ m V V ^ + KV v dpB ^ + ^{iD^ - m^lf) + C h > , (3) 

where my is the resulting mass of the U(l)' vector boson, and Ch> is the Lagrangian for the 
Higgs' particles including their self- interaction and interactions with V^. 

There are four important parameters in the model, the WIMP and vector boson masses 
m^ and my, the mixing parameter k and new coupling constant e'. The most important 
quantity for WIMP physics is arguably the annihilation cross section into the SM states. To 
this end it is easy to identify the two primary mechanisms responsible for annihilation (see 
Fig. [2D: 

(A) ip + -0 — > virtual V — > virtual 7, Z — > SM states. 

(B) xp + xjj — > on-shell V + V, with subsequent decay to SM states. 

Process (B) is open only in the kinematic regime > my while process (A) can occur 
regardless of the relation between my and m^. We will analyze the case of < m v first. 

2.1.1 The characteristic WIMP regime 

The annihilation cross section in this case is given by the diagram in Figure EJ^A). Although 
we use the full result for numerical analysis, to simplify the presentation its helpful to quote 
the annihilation cross section in the limit m 2 z ,mf,m1 m^. Since 2m^, then provides the 
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Figure 2: WIMP annihilation for: (A) < my on the left; and (B) > my on the right - the secluded 
regime in which the annihilation may proceed via two metastable on-shell Vs, which ultimately decay to 
SM states. 



energy scale for the problem, in this limit one may substitute d^B^ v by the total hypercharge 
current and neglect the influence of SM threshold effects. For small mixing, characterized 
by (3 <C 1 where 



Ke' 



\ecosv w J 

the resulting annihilation cross section for nonrelativistic WIMPs takes the following form, 

(a ann ,) m ,» msM - 1.3 pbn x (3 (^^) x [ jT*^ > ( 5 ) 

proceeding in the / = channel with an obvious pole at = my/2, in the vicinity of 
which a more accurate treatment of the thermal average is required. The result depends 
on the mixing parameter j3 and the sum of squares of the hypercharges for the SM fields, 

Sfermions Y f + \ S bosons Y b = 10 + °- 25 - Note that in the opposite limit, m b < my < m z , 
the total cross section is instead proportional to the sum of squares of all the electric charges 
of SM fermions with the exception of the t-quark. 

This cross-section needs to be compared with the constraint on the dark matter energy 
density provided by recent cosmological observations: 

2 x - 1QS W r /) < nDMh 2 ^ 0A (6) 

^JTj) x GeV x M Pl (av) ~ 



where Tf is the freeze-out temperature (it suffices here to take ra^/Tf ~ 20), g>* the effective 
number of degrees of freedom at freeze-out, and the extra factor of two relative to the 
standard formula (see e.g. [16j) is because annihilation can occur only between particles and 
anti-particles. 

In Fig. El we exhibit the abundance constraint on the (3 — plane for a specific choice 
of mediator mass, my = 400 GeV, by saturating the inequality (jSJ). This value of m v 
lies outside the direct reach of LEP or the Tevatron but is certainly within range for the 
LHC. One can clearly see the enhancement of the annihilation cross section in the vicinity 
of the two vector resonance poles, Z and V, where the mixing parameter (3 is allowed to be 
significantly smaller than 1. 

This model is subject to various constraints from direct searches and collider physics. 
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(a) The elastic scattering of galactic WIMPs off nuclei occurs with a characteristic momen- 
tum transfer of order 100 MeV or less, making virtual photon exchange the dominant 
contribution to the scattering amplitude. From an appropriate low-energy effective 
theory viewpoint, the WIMP ip is electrically neutral but exhibits a non-zero charge 
radius given by 

rl = ^*\ = ^ ll2c °: 9w . (7) 



e m v m v 



The contribution of r c to the elastic scattering of WIMPs off nuclei was calculated in 
[T7] . and can easily be rescaled to the equivalent cross section per nucleon, 

47r r?2 2 4 ( m Am D \ 2 _AlT 242 fZ 



°"el = — Z a T c ■ (Tnucleon = — « r m — 

9 \m A + m D J 9 p \A 

where m A and m p are nuclear and nucleon masses, and the Z/A ratio should be specified 
for the relevant experimental setup. We plot the corresponding experimental limit 
recently obtained by the XENON collaboration [3] in Fig. [31 which clearly rules out a 
significant portion of the parameter space away from the = My/2 resonance. This 
is not surprising, as the model is in many ways similar to the original "heavy Dirac 
neutrino" of Ref. [2J, which is known to be essentially excluded by direct searches. 

(b) Other particle physics constraints on this model are also highly dependent on my- For 
my 3> mz and E, where E is the energy accessible in the collision, one can approximate 
K-exchange between SM particles as an effective current- current interaction, 

_ 4W y y 

LeR ~ m v cos 6 2 w » J »- [ } 

We can then constrain the coefficient using limits on the corresponding effective four- 
lepton operator from LEP2 and Tevatron searches for "compositeness" [HI [19]. In 
terms of the conventionally normalized coefficient, Att/A^, no deviation is observed from 
the SM cross section up to A c ~10-15 TeV. Consequently, we arrive at the following 
collider constraint, 

A c <10TeV (3 < 0.3 x — x ( ^ Y . (10) 

1 a V500 GeV/ V ; 

In order to plot this constraint in Fig. [31 we choose two representative values for the 
coupling constant e', one defining the perturbative regime, a' < 1, and a more realistic 
line for a' = a. 

(c) For m^ < mz/2, there is also an extra contribution to the invisible width of Z, namely 

z -> 

a(3 mi r , 2m\ 



Tz^ = ^ 7 2 Z 2 , 2 (l + —?)^m%- 4mJ . (11) 
3 (my — m z z ) 2 \ m z z 1 v 



Requiring this width to be less than 4 MeV [20J results in an additional constraint 
plotted in Fig. [31 
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Figure 3: (A) On the left, the parameter space of the /7(l)'-mediated model. The elastic cross section 
constraint excludes a large portion of the — (3 parameter space. The upper-left corner is also excluded 
due to the invisible decay of Z into WIMPs. The two horizontal lines correspond to the collider limits on 
four-fermion interactions with a' = 1 and a' — a. (B) On the right, the parameter space of the singlet-Higgs 
mediated model. There are no collider constraints, and the <j e \ constraint is much weaker. 

As is evident from the figure, the existing constraints already rule out WIMP masses up 
to 100 GeV with the exception of a small region in the vicinity of m^ = my /2 where there 
is resonant enhancement in the annihilation cross section. Future LHC experiments, and 
the next generation of dark matter experiments, will provide much deeper probes into the 
parameter space of this model. 



2.1.2 The secluded WIMP regime 

The situation changes drastically if annihilation process (B) is kinematically allowed. The 
cross section for WIMP annihilation into pairs of (unstable) V bosons is then given by 



which together with Q implies that in the limit /3 < 1 the correct dark matter abundance 
is achieved if 

1/4 

3 ( m ip 



a ' x i _ ^ ~5x 10~ 3 x '7 . (13) 
^ ml J V500 GeV/ V ; 

This constraint is easily satisfied for a rather natural range for and a'. Crucially, since 
large mixing is no longer required in this kinematic regime to ensure the correct annihilation 
cross-section, k can be taken very small indeed. The only constraints one has to impose are 
that the decay of V (and also h!) occur before the start of nucleosynthesis. By choosing 
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m' h > my/2, only the decays of V are sensitive to the mixing: 

A tighter constraint would follow from requiring V decays to remain in thermal equilibrium, 
which would also ensure the initial thermal and chemical equilibirium for WIMPs, used in 
the derivation of the abundance formula ([H]), 

IV > Hubble Rate [T ~ 0.05m,/,] =► k 2 {j^y) * ^ {^kvY " (15) 

Although considerably tighter than (fT4l) . the constraint ffl5|) does not change the main con- 
clusion: in the limit m me( jiator < wwimp, the WIMP sector can be secluded and neither 
collider nor underground searches impose any significant restrictions on the model. More- 
over, the constraint (|15p can be relaxed to (|14p if some new UV physics ensures proper 
thermal contact between the dark matter and SM sectors at higher temperatures resulting 
in Tdm ~ O(0. 1) x Tsm & t the time of dark matter annihilation. For example, a dimension- 
8 operator of the form V 2 (F^ V ,) 2 generated at a rather high energy scale would suit this 
purpose. 

It is also worth noting that one could completely sever the remaining link to the SM, 
and set the mixing parameter k to zero, if the U(l)' gauge symmetry is not broken at all 
or broken only at very low energy scales. In this case, there will be an extra component 
to the energy density of the universe, namely the "dark radiation" associated with massless 
V^-bosons. However, this energy density associated with V may be much smaller than that 
of the SM photons, because the decoupling of the dark sector may have occurred at a very 
early epoch, after which the photon temperature was effectively increased several-fold by 
input from the decay of numerous SM degrees of freedom. Although at first sight this model 
looks like a completely decoupled dark sector, the energy density in V could in principle be 
detected by highly sensitive next-generation CMB anisotropy probes at very small angular 
scales. 



2.2 Singlet scalar mediator 

An alternative class of models uses scalar mediators, discussed for example in several recent 
publications [TO] . The simplest example consists of two additional fields, the singlet WIMP 
fermion ip and a singlet mediator 0, 

£ W IMP+mediator = ^(^(j)) 2 - -m|0 2 + ^ {id ^ ^ - 77ty - A</,0)^ (16) 

- \ x vcj>{rfH- V ^)-\ 2 <?{rfH- V ^)-V{4>). 

Here H is the SM Higgs doublet and v the corresponding vacuuum expectation value, intro- 
duced in (1171) for convenience. By a redefinition of the field, we can always set (0) = 0. 
Should one pursue the minimalist approach, itself can be a dark matter candidate, which 
would then fix the value of A 2 . However, to demonstrate our point as simply as possible, 
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we take the limit A2 — > 0. After electroweak symmetry breaking, the relevant dark matter 
Lagrangian then takes the following form, 

AviMP+mediator = ~^(d^(f>) 2 + ^{18^^ - - \^4>)^ - \lV 2 (j)h H , (17) 

where all interactions between SM fields and WIMPs are mediated by Higgs-singlet mixing. 

The analysis of this model follows similar lines to the example above, so we will be rather 
brief. We first consider the case > and choose = 400 GeV to enable a clear 
comparison with the previous model. Working to lowest order in the mixing parameter, we 
have computed the annihilation cross section taking into account ti, ZZ, WW, hh and bb in 
the final state, following e.g. Ref. [9]. The cross section has a leading p-wave contribution, 
which provides an additional suppression to (<7f}E Choosing the Higgs mass to be rrih = 120 
GeV, and denoting the mixing parameter by j3h, 

A - (18) 



we can use the abundance constraint ([5J to plot the dependence of j3h on the WIMP mass in 
Fig. 3(B). This plot also includes the experimental constraint on the nucleon-WIMP elastic 
scattering cross section induced by h — <p scalar exchange. We take the SM value of the 
Higgs-nucleon coupling to be ghNN — 300 MeV. There are no significant collider constraints 
on this WIMP model, as all WIMP production has to occur via a real or virtual Higgs or 
Higgs-like bosons, and this is highly suppressed. Although quite constraining for low WIMP 
masses, in this scenario the direct search experiments cannot probe above 50 GeV. 

Again, the situation changes significantly if < m^, which allows for the WIMP sector 
to be secluded. The kinematically allowed pair-annihilation of WIMPs into two 0-scalars 
then simply imposes a constraint on a combination of A^, and m^, but removes any 
constraints on Ai. The decay of these scalars before BBN imposes a very relaxed requirement 
of Ai > 10 -21 . In practice, even taking (3h ~ 10~ 5 would eliminate any chances for direct 
search experiments and/or colliders to probe the WIMP sector in such a model. 



2.3 Right-handed neutrino mediator 

In previous subsections the choice of metastable mediator, although relatively simple, was 
nonetheless exotic. Neither metastable vector or scalar particles are required by any known 
SM physics. There is, however, the distinct possibility of promoting right-handed neutrinos 
Nr, arguably the best motivated extension of the SM field content, to the role of metastable 
mediators with a dark matter sector. Indeed, if the right-handed neutrinos have Majorana 
masses of order the electroweak scale, their decay width into SM states such as left-handed 
neutrinos and the Higgs will be proportional to the square of the Yukawa coupling, i.e. 
in practice to the light neutrino masses. As is clear from the discussion of the previous 
subsections, Tn r ~ m u ~ 0.1 eV is in the right range for Nr to play the role of a metastable 

1 The presence of a pseudoscalar coupling, (jnjji^ijj , would open up the s-channcl. 
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mediator. To complete this model, we choose the secluded sector to have one additional 
fermion N' and a boson S: 



AviMP +med iato r = \{d,Sf -^S 2 + N'i$N> -^fN' T N> + N^Nn-^NlNn 

- ^S 2 H^H - [Y V LHN R - Y N ,SN T R N' + (h.c.)] . (19) 

This Lagrangian possesses a Z 2 symmetry which allows only even numbers of N' or S at 
each vertex. Alternatively, one may complexify S and N' and introduce a new global charge 
that would ensure the same property. 

A viable secluded dark matter model results from the choice mg, ttln' > rriN R . If 
rris > rriN', then S'-mediated annihilation N' + N' — > Nr + Nr will ensure the right WIMP 
abundance of N' given an appropriate choice of Yn>- Notice that a scalar coupling A as 
small as 10~ 8 is sufficient to keep the dark sector in thermal/chemical equilibrium prior to 
freezeout [9]. For rris < tun' the roles are reversed, and A'-mediated annihilation allows for 
S to be dark matter. In either case, the dark matter candidate, either N' or S, is secluded 
from direct observational probes. 



2.4 Strong non-Abelian interactions in the secluded sector 

Models with a non-Abelian gauge group G' in a hidden sector can also be secluded and may 
be of interest in the context of dark matter models with strong self-interactions [15] if the 
scale is relatively light. The gauge bosons of G' cannot couple directly to the SM because of 
gauge invariance. Therefore, in this case the mediators must be charged under both gauge 
groups. The simplest example of this kind is given by 

AviMP+mediator = ~\{G^? + E K^pP^ ""»/)/ + Z) ~ ™*M ( 20 ) 

where is the non-Abelian field strength in the hidden sector, while ip and / are fermions, 
with ip charged only under G', while the fermions / play the role of mediators and are charged 
under both G 1 and the SM gauge group. This field content will necessarily have to satisfy 
anomaly cancelation constraints. We further assume that the confinement scale of G' is 
comparable to or larger than the electroweak scale, A' > v. 

The Lagrangian (TSUI) allows the construction of a secluded WIMP model, in which both 
mediators and dark matter are composite. If all masses in the /-sector are large, these 
fields can be integrated out leading to non-renormalizable interactions between the two 
sectors of the form ^{G a ) 2 (F}^f) 2 . The phenomenology of such terms, in connection 

with metastable dark matter composed of hidden-sector glueballs was considered in Ref. 
|21j . The -^-containing baryons B' and antibaryons are viable dark matter candidates if 
the nonrenormalizable terms leading to their decays are forbidden. The mass of an exotic 
"meson" (M' ~ ipip) would be smaller than the baryon mass, and the annihilation process 
B B' — > M'M' therefore open. Since this cross section is not expected to have additional 
suppression other than that provided by the mass scale of the exotic baryons, the abundance 



9 



constraint would require that the dark matter mass be above 10 TeV. The mesons M' as 
well as the glueballs of the hidden group would be unstable with respect to decay into the 
SM fields, with total widths controlled by the combination A /9 /mj [21] • For large rrif, i.e. 
rrif ^> A', this width can be exceptionally small, e.g. of order the Hubble scale during 
BB freeze-out. Explicit models with exotic baryons as dark matter have previously been 
constructed within the framework of gauge-mediated supersymmetry breaking [22, 23J. 

The opposite limit, rrif <C A', would provide a regime where the WIMP sector could in 
principle be probed because the mesons built from / would then be of electroweak scale, and 
will interact both with the SM and with WIMPs in the form of exotic baryons. Although 
the direct detection signal signal might again be rather low, the exotic /-containing mesons 
could conceivably be produced in the next generation of colliders [25] . 

3 MeV-scale dark matter and mediators 

The choice of a relatively light mediator is perhaps the only viable possibility for WIMP 
masses to lie well below the Lee- Weinberg window, and close to the MeV-scale [23] . This 
mass range has some interesting consequences, including the speculative possibility of linking 
the unexpectedly strong and uniform emission of 511 keV photons from galactic center [7J 
to positrons created by the annihilation 0(MeV)-scale dark matter [5J. To date, much of 
the model-building in this direction has concentrated on utilizing an additional U(l) gauge 
group, under which both the Standard Model and the dark sector are (disproportionately) 
charged: a small charge for the Standard Model fermions, and a larger charge for dark matter 

The most natural anomaly-free quantum number to gauge is B — L [27] , in which case the 
coupling of the additional U(l) gauge bosons to charged fermions is necessarily vector-like. 
The absence of an axial vector current allows this scenario to escape strong constraints from 
atomic parity violation, and from flavor-changing decays of K and B mesons [26]. However, 
there is a price to pay as the coupling to neutrinos creates a problem with the energetics 
of supernovae. During the explosion, the MeV-scale WIMP is thermalized and coupled 
to neutrinos too strongly, supressing the energy of the emitted neutrinos and making the 
observed SN1987A signal highly unlikely [28J. There are two ways to escape this problem: 
taking the WIMP mass in excess of 10 MeV, or forbidding couplings to neutrinos [28J. 
Unfortunately, the first option doesn't work, because the shape of the 511 keV line |29j, 
along with the 7-ray spectrum in the MeV region |30j, do not allow the mass of annihilating 
particles to be in excess of 3-5 MeV. The second option, i.e. no coupling to neutrinos, 
requires abandoning the initial assumption of gauging B — L. 

In this context, it is easy to see that the secluded models of sections 2.1 and 2.2 with 
vector and scalar mediators can solve certain model building problems for MeV-scale dark 
matter. Indeed, neither of these mediation mechanisms lead to any significant coupling to 
neutrinos. In particular, for the U(l)' mediator the if) — u scattering amplitude is necessarily 
suppressed by the Z-boson mass, and since the mixing parameter k is smaller than 1, the 
ip — v scattering cross-section will be below the typical weak-scale value. This motivates 
closer inspection of these two secluded models in the MeV mass range, in order to determine 
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if they do indeed represent viable MeV-scale WIMP scenarios. In what follows we will make 
use of two relations that generalize results of [2H1 EI], namely for the abundance. 



(2-4) 




(21) 



and the 511 keV flux. 
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where mx is the mass of the dark matter candidate X (with 0.5 MeV< nix Izj 3-5 MeV), 
Qx/^dm < 1 is the contribution of X to the total dark matter energy density, and the 
ratio ^li^x/^sii.totai is the contribution of X-annihilation to the total 511 keV 7-ray flux 
from the galactic center. In (|2"2"|) . iV e + is the positron multiplicity per annihilation event. 
Since the neutrino couplings are negligbly small, and direct annihilation to photons is also 
suppressed [32J, in these models e + e~ is the dominant annihilation mode, and N e + = 1 
for mwiMP < "^mediator, and N e + = 2 otherwise. The subscript on the annihilation cross 
section, (av) c or (crv) g , denotes the type of averaging: (-• -) c implies thermal averaging for 
cosmological freezout; while (• • • ) g refers to the average over the galactic velocity distribution 
for the A-particles. We should caution the reader that the relation ( 1221) is only valid to within 
one to two orders of magnitude, due to the uncertainty in the dark matter number density 
in the galactic core. 

(a) U(l) / -mediator, mx > mv- Irrespective of whether the dark matter is a fermion 
or a scalar, annihilation to two V bosons proceeds in the s-channel, and (av) c ~ {cv) g . 
This immediately rules out MeV-scale X particles as the dominant component of dark 
matter, since Qx/^dai — 1 will ensure that the 511 keV 7 ray flux is significantly 
overproduced, ^si^x/^ii^otai ~ 10 4 . Alternatively, one can tune the coupling a' to 
be much smaller than the SM gauge couplings, a' ~ 2 x 10 -3 x a, and satisfy the 
511 keV flux constraint with Qx/^dm ~ 2 x 10~ 5 . This scenario is reminiscent of the 
decaying sterile neutrino [33], where a very small contribution of sterile neutrinos to 
the total dark matter budget can provide the requisite flux. The dominant component 
of dark matter should of course come from other sources, which renders this model 
incomplete. 

(b) U(l) / -mediator, mx < my: In this regime, the choice of X as a scalar charged 
under U{1)' is preferred [24], because of the p-wave annihilation that leads to (<jv) g ~ 
10~ 5 x (av) c . It then appears entirely possible to satisfy both (I2TT) and (|2"2"|) with 
mx ~ O(MeV), Qx/Qdm — 1, and a mixing parameter, 



This value for the mixing parameter would be natural for example if induced radiatively 
at a high scale by a state charged both under the SM U(l) and the extra U(l)'. As 
mentioned before, this model does not pose any problems with respect to the supernova 




(23) 
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signal, and does not presuppose any hierarchy of gauge couplings as a' can be taken 
of order a. Therefore, this model appears the most natural candidate for MeV-scale 
secluded dark matter, having the chance to explain the 511 keV line from the galactic 
center. 

(c) ^-mediator, mx > fn<t>- In this scenario, it is advantageous to have a fermionic 
dark matter candidate ip with scalar (rather than pseudoscalar) couplings to 0. The 
annihilation ipip — > <ft<j) proceeeds in the p-wave and can always be tuned to the required 
level with a typical choice ~ 10 -6 . Since ~ few MeV, this value of the Yukawa 
coupling is natural. The subsequent decay of due to mixing with the Higgs is highly 
suppressed by the electron Yukawa coupling, 

x (j^y x ^ ^ (^) 2 >io- (24) 

m h J WW 8tt V K / 

The naturalness requirement for the 0-mass would impose a significant constraint here. 
If we consider the contribution from Higgs mixing in (fT7|) . Xxv/rrih < m^/v, this clearly 
favors a long 0-lifetime (~ 1 sec) and a small mixing parameter. Even then, one must 
ensure that the "missing energy" decay K + — > 7r + + is within the allowed range. At 
the quark level, the amplitude for the process is given by a Higgs penguin (see, e.g. 

mr- 

£eff = [y^r ) 6A7r2m> d L 8 R 4> + (h.c), (25) 
leading to the (non-SM) missing energy decay, 



AiuV f3m? t V td V*\ 2 m\ 



K—>ir+(h— mediator — I oil -,r>99 I n a •}' 



(26) 



Requiring that this width not exceed the observed missing energy decay branching 
ratio Br = 1.5]tJJ x 10 -10 [35] associated with the SM process K + — > nuu, results in 
the following constraint on — h mixing: 

^) 2 <2xl0-. (27) 
m h / 

This cuts out a significant part of the parameter space, but together with (1241) still 
leaves a relatively narrow interval for the mixing parameter, 10~ 7 — 10 -8 , where the 
model survives all constraints (although not without a modest amount of fine-tuning 
of the mediator mass) and thus can be the dominant dark matter component while still 
accommodating the positron signal through a combination of annihilation and decay. 

The constraints remain essentially the same for a pseudoscalar coupling of (f) to the 
fermion ip, if the Higgs sector in SM is assumed to be minimal, in which case the 
mixing constant Ai is CP-violating. The additional processes: s-wave annihilation 
if)if> — > e + e~ through a virtual 0, and also ipip — > (jxficj) if kinematically allowed, are too 
weak in comparison with the p-wave annihilation ipip — > cfxp to affect the constraints 
discussed above. In principle, with an extended Higgs sector, could also mix in a 
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CP-conserving way with the physical CP-odd Higgs scalar (s) and in addition could 
have an enhanced coupling to electrons, relative to the top quark, thus significantly 
relaxing the constraints on the parameter space in comparison with Eqs. (I24p and (I27p . 

(d) 0-mediator, mx < m^'- In this case the annihilation cross section is suppressed 
by (m e /v) 2 ~ 4 x 10~ 12 and the relation ( 12~TT) would require a mixing parameter, fih, 
of order one. Such a choice would involve a gross violation of naturalness, and would 
also lead to an unacceptably large missing energy signal in K and B decays [34j. We 
therefore conclude that this option is not viable. 

4 Concluding Remarks 

Secluded WIMP dark matter appears to be a generic possibility, as rather minimal model- 
building choices lead to viable WIMPs interacting with metastable mediators. These me- 
diators could be either elementary or composite, and we have constructed explicit models 
with scalar, vector and fermion mediation to the Higgs, hypercharge gauge boson, and light 
neutrino sectors respectively of the SM. In the latter case there is a natural choice for the 
mediator, namely a right-handed electroweak-scale neutrino. 

Despite existing as a thermal relic, with weak-scale annihilation, none of the secluded 
WIMP models constructed here lead to appreciable signals in underground detectors or 
register as a missing energy channel in collider experiments. Nonetheless, these models 
are subject to indirect constraints related to 7-rays caused by WIMP annihilation, e.g in 
the galactic center. Therefore, one of the main conclusions of this paper is the complete 
complementarity of direct and indirect efforts for detecting WIMP dark matter. 

Two of the WIMP models constructed here allow a rescaling down to the MeV mass range, 
motivated by the intriguing connection to the galactic 511 keV line. The secluded models 
considered here do not couple dark matter to neutrinos, and therefore avoid problems with 
the suppression of SN1987A signal. We have shown that models with an additional U(l)' 
and kinetic mixing to the hypercharge gauge boson circumvent some problems of scenarios 
where SM fields carry an additional gauge quantum number, and appear to be free from 
unnaturally small parameters. The main limitations come instead from the dark matter 
energy density and from indirect constraints related to the galactic 7-ray spectrum in the 
MeV region. 



Note added — As this paper was being finalized, we became aware of a recent preprint [36] 

that also deals with U(l)' models of MeV-scale dark matter with kinetic mixing, and thus 
has some overlap with the discussion in Sect. 3(a,b). 

Acknowledgements 

The work of MP and AR was supported in part by NSERC, Canada. Research at the 
Perimeter Institute is also supported in part by the Government of Canada through NSERC 
and by the Province of Ontario through MEDT. The work of MV is supported in part by 
the DOE grant DE-FG02-94ER-40823. 




13 



References 



[i] 

[2] 

[3] 
[4] 
[5] 
[6] 
[7 



[9 



G. Fuller et al., arXiv:nucl-ex/0702031 



B. W. Lee and S. Weinberg, Phys. Rev. Lett. 39, 165 (1977); M. I. Vysotsky, A. D. Dol- 
gov and Y. B. Zeldovich, JETP Lett. 26, 188 (1977) [Pisma Zh. Eksp. Teor. Fiz. 26, 
200 (1977)]. 

D. S. Akerib et al. [C DMS Collaboration], Phys. Rev. Lett. 96, 011302 (2006) 
|arXiv:astro-ph/0509259] . 

J. Angle et al. [XENON Collaboration], Dark Matter Experiment at the Gran Sasso 
larXiv:0706.0039l [astro-ph]. 

E. Moulin [H.E.S.S. Collaboration], larXiv:0710.2493 [astro-ph]; A. Malinin [AMS Col- 
laboration], AIP Conf. Proc. 903, 638 (2007). 

C. Boehm, D. Hooper, J. Silk, M. Casse and J. Paul, Phys. Rev. Lett. 92, 101301 (2004) 
|arXiv:astro-ph/0309686] ; 

J. Knodlseder et al, Astron. Astrophys. 411, L457 (2003) [arXiv:astro-ph/0309442| ; 
P. Jean et a/., Astron. Astrophys. 407, L55 (2003 ) [arXiv:astro-ph/03 09484] ; J. Kn- 
odlseder et al, Astron. Astrophys. 441, 513 (2005) |arXiv:astro-ph/0506026|. 

V. Silveira a nd A. Zee, Phys. Lett . B 161, 136 (1985); J. McDonald, Phys. Rev. D 50, 
3637 (1994) |arXiv:hep-ph/0702143] ; 

C. P. Burgess, M. Po spelov and T. ter Veldhuis, Nucl. Phys. B 619, 709 (2001) 
|arXiv:hep-ph/0011335| . 



[10] H. Davoudiasl, R. Kita no, T. Li and H. Murayama, Phys. Lett. B 609, 117 (2005) 
[arXiv:hep-ph/0405097|; C. Bird, R. Kow alewski and M. Pospelov, Mod. Phys. Lett. 
A 21, 457 (2006) |arXiv:hep-ph/0601090]; Y. G. Kim and K. Y. Lee, Phys. Rev. D 
75, 115012 (2007) |arXiv:hep-ph/0611069]; D. P. Finkbeiner and N. Weiner, Phys. 
Rev. D 76, 083519 (2007) [arXiv:astro-ph/070258^ |; F. D'Eramo, Phys. Rev. D 76, 
083522 (2007) [arXiv: 0705 .44931 [hep-ph]]; V. Barger, P. Langacker, M. McCaskey, 
M. J. Ramsey- Musolf and G. Shaughnessy, larXiv:0706.4311l [hep-ph]; H. Sung Cheon, 
S. K. Kang and C. S. Kim, i arXiv:0710.24T6l [hep-ph]. 

[11] S. Dodelson and L. M. Widrow, Phys. Rev. Lett. 72, 17 (1994) |arXiv:hep-ph/9303287| ; 
A. D. Dolgov and S. H. Hansen, Astropart. Phys. 16, 339 (2002) 
|arXiv:hep^ph /0009083|. 

[12] M. Cirelli, N. Forne ngo and A. Strumia, Nucl. Phys. B 753, 178 (2006) 
[arXiv:hep-ph/0512090|; R. Barbieri , L. J. Hall and V. S. Rychkov, Phys. Rev. D 74, 
015007 (2006) |arXiv:hep-ph/060318~8 |; L. Lopez Honorez, E. Ne zri, J. F. Oliver and 
M. H. G. Tytgat, JCAP 0702, 028 (2007) |arXiv:hep-ph/0612275| ; M. Cirelli, A. Stru- 
mia and M. Tamburini, Nucl. Phys. B 787, 152 (2007) [ arXiv:0706.407ll [hep-ph]]. 



14 



arXiv:hep-p h/0302215|. 

P. N. Spergel and P. J. Steinhardt, Phys. Rev. Lett. 84, 3760 (2000) 
|arXiv:astro-ph/9909386| . 

E. W. Kolb and M. S. Turner, The Early Universe, (Westview Press, 1990). For a recent 
discussion of generic WIMP physics see, e.g. M. Drees, H. Iminniyaz and M. Kakizaki, 
larXiv:0704. 1590 [hep-ph]. 



M. Pospelov and T. ter Veldhuis, Phys. Lett. B 480, 181 (2000) |arXiv:hep-ph/0003010|. 



[13] B. de Carlos and J. R. Espinosa, Phys. Lett. B 407, 12 (1997) | arXiv:hep-ph/9705315| ; 
V. Barger, C. Kao, P. Langacker and H. S. Lee, Phys. Lett. B 600, 104 (2004) 
|arXiv:hep-ph/0408120|; D D. F eldman, B. Kors and P. Nath, Phys. Rev. D 75, 023503 
(2007) [arXiv:hep-ph/0610133l; D. Fe ldman, Z. Liu and P. Nath, Phys. Rev. D 75, 
115001 (2007) larXiv:hep-ph/0702l23 |; K. Cheung and T. C. Yuan, JHEP 0703, 120 
(2007) jarXiv:hep-ph/0701107]; W. F~ Chang , J. N. Ng and J. M. S. Wu, Phys. Rev. 
D 75, 115016 (2007) |arXiv:hep-ph/0701254] ; G. Belanger, A. Pukhov and G. Ser- 
vant, larXiv:0706.0526l [hep-ph]; J. H. Huh, J. E. Kim, J. C. Park and S. C. Park, 
larXiv:0711.3"528l [astro-ph]. 

[14] J. R. Ellis, D. V. Nanopoulos and S. Sarkar, Nucl. Phys. B 259, 175 (1985); J. R. Ellis, 
G. B. Gelmini, J. L. Lopez, D. V. Nanopoulos and S. Sarkar, Nucl. Phys. B 373, 399 
(1992); J. L. Feng, A. R ajaraman and F. Takayama, Phys. Rev. Lett. 91, 011302 (2003) 

[15 
[16 

[17 
[18 
[19 
[20 
[21 
[22 

[23 

[24 

[25 
[26 

[27 
[28 



LEP EWWG, ff Subgroup, |http : //lepewwg . web . cern . ch/L EPEWWG/ lep2 



M. Titov [DO Collaboration], PoS HEP2005, 314 (2006) |arXiv:hep-e^7051 2006 ] . 

PDG Review of Particle Properties, W.-M. Yao et al., Journal of Physics G 33, 1 (2006). 

A. E. Faraggi and M. Pospelov, Astropart. Phys. 16, 451 (2002) [arXiv:hep-ph/0008223| . 

S. Dimopoulos, G. F . Giudice and A. Pomarol, Phys. Lett. B 389, 37 (1996) 
|arXiv:hep-ph/9607225| . 

K. Hamaguchi, S. Shirai and T. T. Yanagida, Phys. Lett. B 654, 110 (2007) 
jarXiv:0707.2~463l [hep-ph]]. 

C. Boehm and P. Fayet, Nucl. Phys. B 683, 219 (2004) |a rXiv:hep-ph/0305261| ; 
C. Boehm, P. Fayet and J. Silk, Phys. Rev. D 69, 101302 (2004) |arXiv:hep-ph/0311143|. 



M. J. Strassler and K. M. Zurek, Phys. Lett. B 651, 374 (2007) |arXiv:hep-p h/0604261|. 

P. Fayet, Phys. Rev. D 74, 054034 (2006) |arXiv: hep-ph/ 06 07318] ; P. Fayet, Phys. Rev. 
D 75, 115017 (2007) [arXiv:hep-ph/0702176] . 

P. Fayet, Nucl. Phys. B 347, 743 (1990). 

P. Fayet, D. Hoop er and G. Sigl, Phys. Rev. Lett. 96, 211302 (2006) 
|arXiv:hep-ph/0602169| . 

15 



[29] P. Jean et al, Astron. Astrophys. 445, 579 (2006) larXiv:astro-pli70 509298 1 . 

[30] J. F. Beacom and H. Yuksel, Phys. Rev. Lett. 97, 071102 (2006) 
[arXiv: astro-ph /05 1 24 1 1| . 

[31] Y. Ascasibar, P. Jean, C. Boehm and J. Knoedlseder, Mon. Not. Roy. Astron. Soc. 368, 
1695 (2006) [arXiv:astro-ph/050714 2| . 

[32] M. Pospelov and A. Ritz, Phys. Lett. B 651, 208 (2007) ||arXiv:hep-ph/0703128| . 

[33] C. Picciotto and M. Pospelov, Phys. Lett. B 605, 15 (2005) | iarXiv:hep-ph/0402178] . 

[34] C. Bir d, P. Jackson, R. Kow alewski and M. Pospelov, Phys. Rev. Lett. 93, 201803 
(2004) [arXiv:hep-ph/0401195] . 

[35] S. Adler et al. [E787 Collaboration], Phys. Rev. Lett. 88, 041803 (2002). 

[36] J. H. Huh, J. E. Kim, J. C. Park and S. C. Park. larXiv:0711.3528l [astro-ph]. 



16 



